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We report a measurement of the exclusive e'^e" -^ A^A^T cross section as a function of center- 
of-mass energy near the A+A~ threshold. A clear peak with a significance of 8.2cr is observed in 
the AJA^ invariant mass distribution just above threshold. With an assumption of a resonance 
origin for the observed peak, a mass and width of M = (46341* (stat.)^g(sys.)) MeV/c^ and Ftot = 
(92l24(stat.)_2^(sys.)) MeV are determined. The analysis is based on a study of events with initial- 
state-radiation photons in a data sample collected with the Belle detector at the T(4S') resonance 
and nearby continuum with an integrated luminosity of 695 fb~^ at the KEKB asymmetric-energy 
e"'"e~ collider. 
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The discovery of many unexpected charmonium-likc 
states has stimulated renewed interest in charmonium 
physics. Among these new states, the y(4260) [l|, 0j 
y(4360) and y(4660) 0, 01 have quantum numbers 
J^*-^ — 1 and are produced via e~^e~ annihilation. 
Surprisingly, no evidence for open-charm production as- 
sociated with these new states has been observed. More- 
over, the parameters of the conventional charmonium 
1~~ states obtained from fits to the inclusive cross sec- 
tion [^ remain poorly understood theoretically Q . Mea- 
surements of exclusive cross sections for charmed meson 
and baryon pairs in the 4 to 5 GeV energy range are 
needed to help to clarify the situation. 

Initial-state radiation (ISR) provides a powerful tool 
for measuring exclusive e'^e^ cross sections at ^/s smaller 
than the initial e^e^ centcr-of-mass (cm.) energy 
(^c.m.) at i?- factories. ISR allows one to obtain cross 
sections over a broad energy range, while the high lu- 
minosity of the i?-factories compensates for the sup- 
pression associated with the emission of a hard pho- 
ton. The first measurements of the exclusive cross sec- 
tions for e+e" -^ d'^*)±d*^ for ^ near the D^*'>^D*^ 
thresholds were performed by Belle [Tj. Subsequently, 
BaBar [8| and Belle [9(] presented exclusive e~^e~ — > DD 
production measurements via ISR. Recently, Belle [10| 
reported a measurement of the exclusive cross section 
for e~^e^ -^ D^D~Tr~^ [ll| and the first observation 
of V'(4415) -^ DD^{2A60) decay. These measured final 
states almost saturate the total cross section for hadron 
production in e"'"e~ annihilation in the ^/s region up to 
~ 4.3 GeV. The thresholds for charm baryon-antibaryon 
pair production lie in the energy range above 4.5 GeV, 



where experimental data are limited [l^l or unavailable. 
In this Letter we report the first measurement of the 



exclusive cross section for the process e"^ 



A+A: 



via ISR and the first observation of a resonant-like struc- 
ture at threshold. The data sample corresponds to an 
integrated luminosity of 695 fb~ collected with the Belle 
detector [l3| at the T(45') resonance and nearby contin- 



uum at the KEKB asymmetric-energy e+e" collider [14|. 
The selection of e+e~ — > A+A~7isr signal events using 
full reconstruction of both the A+ and A~ baryons suf- 
fers from the low Ac reconstruction efficiency and small 
branching fractions for decays to accessible final states. 
Therefore, in order to achieve higher efficiency we require 
full reconstruction of only one of the Ac baryons and the 
7ISR photon. In this case the spectrum of masses recoil- 
ing against the A+71SR system. 



M,ec(A+7iSR) = 



(^c.m. — ^!+ ) 

Ad7lSR 



|2_„*2 



P 



A+7isr' 



(1) 



peaks at the Ac mass. Here E* 



and p* 



7.1 cxii»a iJ . .^ are 

A?7isR A;r7isR 

the center-of-mass energy and momentum, respectively, 
of the A+71SR combination. The Mi-oc(A+7isr) peak is 
broad {<jm^^^ ^ 250 McV/c^ according to a Monte Carlo 
(MC) simulation) and asymmetric due to the photon en- 
ergy resolution and higher-order ISR processes {i.e. more 
than one 7isr in the event). This makes the distinction 
between A+Aj, A+Ajtt" and A+Ajtttt final states diffi- 
cult. 

For the measurement of the exclusive cross section 
for e+e^ — > A^Aj, we determine the mass recoiling 
against the 71SR photon (Mrec(7iSR))5 which is equivalent 
to M(A+A~) in the absence of higher-order QED pro- 



cesses. To improve the Miec(7isR) resolution (expected 
to be ~ lOOMeV/c^), we apply a refit that constrains 
M].ec(A+7isR) to the nominal A~ mass. In this way 
we use the well measured properties of the fully recon- 
structed A+ to correct the poorly measured energy of the 
7iSR. As a result, the M^+^- resolution is improved sub- 
stantially; it varies from ~ 3 MeV/c^ just above threshold 
to ~ 8MeV/c2 at Mj^+j^- ~ 5.4GeV/c2. 

All charged tracks are required to originate from the 
vicinity of the interaction point (IP); we impose the re- 
quirements dr < 1 cm and jdzj < 4 cm, where dr and \dz\ 
are the impact parameters perpendicular to and along the 
beam direction with respect to the IP. Particle identifica- 
tion requirements are based on dE/dx, aerogel Cherenkov 
and time-of-flight counter information [l5|. Protons and 
charged kaons have typical misidentification probabili- 
ties less than 0.1. No identification requirements are ap- 
plied for pion candidates. K^ (A) candidates are recon- 
structed from tt^t:^ (p?:^) pairs with an invariant mass 
within lOMeV/c^ (~ 3ct) of the K'^ (A) mass. The dis- 
tance between the two pion (proton and pion) tracks at 
the Kg (A) vertex must be less than 1 cm, the transverse 
flight distance from the interaction point is required to 
be greater than 0.1 cm, and the angle between the Kg{A) 
momentum direction and the flight direction in the x — y 
plane should be less than 0.01(0.005) rad. Photons are re- 
constructed in the electromagnetic calorimeter as show- 
ers with energies greater than 50 MeV that are not asso- 
ciated with charged tracks. ISR photon candidates are 
required to have energies greater than 3.5 GeV. Candi- 
date 7r° mesons are formed from pairs of photons. If the 
mass of a 77 pair lies within 15MeV/c^ (^ 3(t) of the 
tt" mass, the pair is fit with a tt'^ mass constraint and 
considered as a 7r° candidate. 

A+ candidates are reconstructed using three decay 
modes: pKg, pK~Tr^ and Att"*". The mass distribution 
of A+ candidates from A+71SR combinations is shown in 
Fig. [T](a). To suppress combinatorial background, we re- 
quire the presence of at least one p in the event from the 
decay of the unreconstructed A~ {p tag). As a result, 
the combinatorial background is suppressed by a factor 
of ~ 10 at the expense of about a 40% reduction in signal 
according to the MC simulation (see Fig.[l](b)). 

A ±10MeV/c^ mass window is used for all A+ can- 
didate decay modes (~ 2.5 a in each case). To improve 
the momentum resolution of Ajj" candidates, final tracks 
are fitted to a common vertex with a mass constraint to 
the A+ mass. Only one A+71SR combination per event 
is accepted; in the case of multiple combinations, which 
occur in 5% of the candidate events, the combination 
with the best x^ for the A+ mass fit is selected. A+ 
mass sidebands selected for the background study are 
four times as large as the signal region. To avoid signal 
over-subtraction, the sidebands are shifted by 20 MeV/c^ 
from the signal region. The sidebands are divided into 
windows of the same width as that for the signal. The A+ 
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FIG. 1; The mass distribution of Aj from AJ71SR combina- 
tions: (a) without a p tag; (b) with a p tag. The Aj signal 
and background regions are indicated by vertical lines. 



candidates from these sidebands are refitted to the cen- 
tral mass value of each window and a single candidate in 
each window per event is selected. 

The distribution of M^cc (A^7isr) with a p tag is shown 
in Fig. [21 The excess around the A~ mass includes the 
e+e~ — > A+A~7isr signal as well as possible reflections 
from the e'^e^ — > A+A~7r'^7isi. and e'^e" -^ A'^ A'^ mryisi: 
processes with an additional tt*^ or tttt, respectively, in 
the final state. The process e+e^ 
could proceed via e+e^ -^ A+S 
and is expected to be strongly suppressed. The pro- 
cess e~^e~ — > A+Aj7r7r7isi- is allowed and is expected to 
proceed via A+ A" (2595), A+A^(2625), A+A-(2765) and 
A+Aj(2880) final states. Each final state would produce 
a broad peak in the Mrec(A+7isR) distribution around 
the corresponding mass value {i.e. m^- , rn/^-(^2595)^ 
n^A-(2625), "^A-(2765) and toa-(2880))- Due to the poor 
Afi.cc(A^7iSR) resolution these peaks overlap and appear 
as a shoulder for masses above ^ 2.5GeV/c^. 

To estimate the contribution from the reflections and 
to optimize the signal region requirement, we fit the 
A'/rcc(A+7isR) distribution with the sum of a signal plus 
a combinatorial and reflection background with normal- 
izations left as free parameters. To describe the combi- 
natorial background, we use A+ sideband data param- 
eterized by a second-order polynomial. We perform a 
simultaneous likehhood fit to the Mrec(A+7isR) signal 
and sideband spectra. The signal and reflection shapes 
of the A+S-, A+A-(2595), A+Aj(2625), A+Aj(2765), 
A+A~(2880) final states are fixed from the MC sim- 
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FIG. 2: The Mrec(A+7isR) distribution with a p tag. The 
solid curve represents the result of the fit described in the text. 
The combinatorial background parameterization is shown by 
the dashed curve. The dashed-dotted curve represents a con- 
tribution of the AjA~ final state while the dotted one is that 
of the K+K' (2595) and the A+ A;: (2625) final states. The dif- 
ference between the solid and dotted lines corresponds to the 
contribution of the A+A;r(2765) and the A+A;:(2880) final 
states. The histogram shows the normalized M .+ sidebands 
contributions. The selected signal window is indicated by the 
vertical lines. 



ulation. All reflection normalizations are floated sep- 
arately in the fit. The goodness of the fit is found 
to be x^/n.d./ = 18.8/22. We define an asymmet- 
ric requirement on A/rcc(A+7isR) of — 250MeV/c^ < 
771^- < 150MeV/c^ to suppress the dominant part of 
the reflection background, as shown in Fig. [2l We flnd 
(386 ± 27(stat.)) signal events in this signal region. The 
contribution of the process e"'"e^ -^ A^A~7r''7isr in the 
signal region is estimated to be less than 18 events at 
the 90% C.L. while that from the e+e^ -^ A+A^7r7r7isi- 
process is estimated to be (7.3± 1.7(stat.)) events. In the 
following study the possible contribution of these back- 
grounds is included in the systematic error. 



The contribution from e^t 



A+A; 



where an 



energetic tt*^ is misidentified as a single 7isr, is found to 
be negligibly small. This is determined from a study of 
e"'"e~ -^ A+A|77r*' events using a similar reconstruction 
technique, but with an energetic tt" replacing the 7isr. 

The Mj^+j^- spectrum for events in the signal region 
is shown in Fig. [3](a). A clear peak is evident near 
the A+Aj threshold. We perform a simultaneous like- 
lihood fit to the M^+^- distributions for the A+ signal 
and sideband regions to fix the combinatorial background 
shapes. The combinatorial background is parameterized 
by piVM-Mthr • e-(P^*^+P^-^^'), where pi, p2 and ps 
are free parameters. The signal function is a sum of a rel- 
ativistic s-wave Breit-Wigner (RBW) function [l6| and a 
threshold function y/M — Mthr with a fioating normaliza- 
tion to take into account a possible non-resonant contri- 
bution. Finally, the sum of the signal resonance and non- 
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FIG. 3: The M^+^- spectrum for the signal region: (a) with 
p tag. The solid curve represents the result of the fit de- 
scribed in the text. The threshold function is shown by the 
dashed curve. The combinatorial background parameteriza- 
tion is shown by the dashed-dotted curve; (b) with proton 
(wrong-sign) tag. Histograms show the normalized contribu- 
tions from At sidebands. 



resonant functions is multiplied by an efficiency function 
that has a linear dependence on M^+y^-, and the differen- 
tial ISR luminosity, described in Ref. [7|]. The fit, shown 
as a solid curve in Fig. [3](a), attributes 142t28(stat.) 
events to the RBW signal. The obtained peak mass 
is M ^ (4634l7(stat.):^8(sys.))MeV/c2 and the total 

width is Ttot = (92l^°(stat.)+2?(sys.))MeV. The fit 
gives 'x^ jn.d.] = 104/77. Here the systematic uncer- 
tainties are obtained by varying the fit range, histogram 
bin size, efficiency function, parameterization of the back- 
ground function and the non-resonant parametrization. 
The systematic error associated with the possible inter- 
ference between the resonance and non-resonant contri- 
butions is estimated from the fit with a coherent sum of 
the RBW and non-resonant amplitudes, which has the 
quality y^ /n.d.f = 103/76 and yields a smaller mass 
(4626MeV/c2) and total width (77MeV). A statistical 
significance for the signal of 8.8(7 is determined from the 
quantity — 2 ln(£o/'Cmax), where £max is the maximum 
likelihood returned by the fit, and Cq is the likelihood 
with the amplitude of the Breit-Wigner function set to 
zero, taking the reduction in the number of degrees of 
freedom into account. The significance including system- 
atics is 8.2(7. We use X(4630) to denote the observed 
structure. 



As a cross check, we present in Fig. [3](b) the A/ 



a+a; 



spectrum for the signal region for wrong-sign tags, i.e. 
requiring a presence of a proton in the event in addition 
to the A+7igR combination. The AIj^+j^- distribution 
from the signal A+ window is in good agreement with 
the normalized contributions from the A+ sidebands. 

The e+e^ -^ A+Aj cross section is extracted from 
the background-subtracted A+Aj mass distribution fol- 
lowing the procedure described in Ref. [7[, taking into 
account the differential ISR luminosity and the efficiency 
function. The resulting e^e~ —> A+A~ exclusive cross 
section is shown in Fig. 5] with statistical uncertainties 
only. Since the bin width is much larger than resolution, 
no correction for resolution is applied. 



TABLE I: Contributions to the systematic error on the cross 
sections, [%]. 



Source 


A+A- 


Background subtraction 


±7 


Cross section calculation 


±3 


Reconstruction 


±5 


Identification 


±3 


Angular distributions 


±4 


Total 


±10 


B(Ac+) 


±26 


B{A- ^ pX) 


±32 
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FIG. 4: The cross section for the exclusive process e'^e 
A+A7. 



The peak cross section for the e+e -^ A+A^ process 
at -E'c.m. = '7ix(463o) is calculated from the amplitude of 
the RBW function in the fit to be a{e+e- -^ X(4630)) x 
S(X(4630) -. A+A-) = (0.47triJ(stat.)+°:[;^(sys.) ± 
0.19(sys.)) nb. Here the first systematic uncertainty is 
obtained by varying the fit range, histogram bin, pa- 
rameterization of the background function, efficiency and 
the possible interference between the resonance and non- 
resonant contributions. The second one comes from the 

= (5.0±1.3)-10~2and 
jig . Using a{e~^e~' — > 



uncertainties in S(A+ -^ pK tt"' 

S(A- -^ pX) = (50 ± 16) • 10- 

X(4630)) = 127r/m2^(4g3o) x (Fee/Ftot) and the X(4630) 

mass value obtained from the fit we calculate Fge/Ftot x 

S(X(4630) ^ A+A-) = (0.68l°:i^(stat)+°-°I(sys.) ± 

0.28(sys.)) • 10-*^. 

The various contributions to the systematic errors for 
the a{e^e^ — > A^A^r) measurements are summarized 
in Tabic HI The systematic errors associated with the 
combinatorial background subtraction arc estimated to 
be 3% due to an uncertainty in the scaling factors for 
the sideband subtractions. It is estimated using fits 
to the Mj^+ distribution with different signal and back- 
ground parameterizations. Reflections from the e'^e~ —> 



A+A; 



"7r°7isr 



and e~^t 



A+A: 



7r7r7isi. processes are es- 



timated conservatively to be smaller than 6% of the sig- 
nal. The uncertainty due to a possible e+e" -^ A+Aj7r° 
contribution is found to be 1%. The systematic error as- 
cribed to the cross section calculation includes a 1.5% er- 
ror on the differential luminosity and 2% error due to the 
MC statistics. Another source of systematic error comes 
from uncertainties in track and photon reconstruction ef- 
ficiencies (1% per track and 1.5% per photon). Another 
contribution comes from the uncertainty in the kaon and 
proton identification efficiency. The systematic uncer- 
tainty due to the unknown helicity angle distribution for 
the A+A" final state is included. For the efficiency cal- 
culation we use a flat helicity distribution and consider 
the extreme cases dN/dcos0 ~ 1 + cos^9 and ~ sin^6' for 
the efficiency uncertainty. 

In summary, we report the first measurements of the 
e+e^ -^ A+A^ exclusive cross section over the center-of- 
mass energy range from the threshold to 5.4 GeV with 
initial-state radiation. We observe a significant near- 
threshold enhancement in the studied cross section. The 
nature of this enhancement remains unclear. In many 
processes including three-body B meson baryqnic de- 
cays, mass peaks are observed near threshold [l7|. How- 
ever, the cross section for e~^e~ — > AA measured via 
ISR by BaBar [l8| has a different pattern: it increases 
sharply at threshold and then decreases gradually with- 
out any peak-like structure. Assuming the observed 
peak to be a resonance, its mass and width are found 
to be M = (4634+?(stat.)+g(sys.))MeV/c2 and Ftot = 
(92l24(stat.)+2°(sys.))MeV, respectively. These values 
are consistent within errors with the mass and width of 
a new 1 charmonium-likc state, the Y(4660), that was 
found in iIj{2S)'ktt decays via ISR [4|. Finally, we cannot 
exclude the possibility that the observed enhancement is 
the 5'^ 5*1 charmonium state that is predicted around the 
observed mass [19[. 
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